Open heavy flavor at RHIC and LHC in a partonic 
transport model* 

Jan Uphoff^, Oliver Fochler^, Zhe Xu^, and Carsten Greiner 

^ Institut fiir Theoretische Physik, Johann Wolfgang Goethe-Universitat 
Frankfurt, Max-von-Laue-Str. 1, D-60438 Frankfurt am Main, Germany 

^ Physics Department, Tsinghua University, Beijing 100084, China 



Heavy quarks are a unique probe to study the medium produced in 
ultra-relativistic heavy ion colhsions. Within the partonic transport model 
Boltzmann approach to multi-parton scatterings (BAMPS) the production, 
energy loss and elliptic flow of heavy quarks are investigated in 3+1 dimen- 
sional simulations of the quark gluon plasma evolution. With only binary 
interactions between heavy quarks and particles from the medium calcu- 
lated within perturbative QCD, the results on elliptic flow and the nuclear 
modification factor are not compatible with experimental data from RHIC 
and LHC. However, if the binary cross section is multiplied with K = 4 
both the elliptic flow and the nuclear modification factor are simultaneously 
described at RHIC and also LHC. Furthermore, preliminary results are pre- 
sented that the implementation of radiative processes leads to a stronger 
suppression which agrees well with the measured nuclear modification fac- 
tor at RHIC without the need of any K factor. 

1. Introduction 

Several evidences [l| show that in ultra-relativistic heavy ion collisions a 
unique state of matter is produced. This medium consists only of deconfined 
quarks and gluons and is, therefore, called quark gluon plasma. Since the 
plasma itself cannot be seen directly, its properties must be revealed by 
measuring particles produced after hadronization of the deconfined medium. 
Due to these indirect and complex measurements, many open questions 
about the properties of the produced medium remain. 

Heavy quarks, that is, in this context, charm and bottom quarks, are 
a perfect probe to tackle many of the open questions. They are well cali- 
brated in a sense that they are produced entirely in the early stage of the 
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heavy ion collision due to their large mass [2] and are also tagged during 
hadronization due to flavor conservation. Charm (bottom) quarks hadronize 
to D (B) mesons, which are, however, difficult to measure. Therefore, at 
RHIC one measures electrons stemming from the meson decays, which can 
also reveal much information about heavy quarks. However, as a disadvan- 
tage it cannot be distinguished between electrons from charm and bottom 
quarks. Nevertheless, at LHC for the first time it is possible to reconstruct 
D mesons and, therefore, receive information only about charm quarks with 
- in principle - no bottom contribution. 

The heavy flavor electron data from RHIC [1, 0] and the heavy flavor 
electron, muon and D meson data from LHC show that the suppression 
of heavy quarks is comparable to that of light quarks. From the theory 
perspective it was thought that radiative processes involving heavy quarks 
are suppressed due to the dead cone effect 0, S], which would imply a 
smaller suppression for massive particles. Elliptic flow V2 measurements of 
particles associated with open heavy flavor also show that heavy quarks 
interact strongly with the other particles of the medium. Whether these 
observations can be explained by collisional or radiative energy loss or other 
effects is currently in debate [§, [13] • 

The paper is organized as the following. After a short introduction to 
our model BAMPS, we compare in Sec.[3]the results on the elliptic flow and 
nuclear modification factor to experimental data from RHIC and LHC. How 
this picture is altered by introducing radiative collisions of heavy quarks 
in addition to binary interactions is shown in Sec. HI followed by a short 
summary. 



2. Partonic transport model BAMPS 

For the simulation of the QGP we use the partonic transport model 
Boltzmann Approach to MultiParton Scatterings (BAMPS) [11], which de- 
scribes the full space-time evolution of the QGP by solving the Boltzmann 
equation for on-shell partons and pQCD interactions. The following pro- 
cesses are implemented for gluons (g) and heavy quarks (Q): gg gg, 
99 999, 999 ^ 99, 99 ^ QQ, QQ ^ 99, 9Q ^ dQ and gQ -\9Q- De- 



tails of the model and the employed cross sections can be found in [3, lid . Ill | . 
In Sec. m we will also discuss some preliminary results including the radiative 
process gQ — > gQg. 



3. Elliptic flow and suppression of heavy quarks 

The elliptic flow V2 and the nuclear modification factor Raa of heavy 
quarks are important observables. Figure [T] compares our results to the 
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Figure 1. Elliptic flow V2 (left) and nuclear modification factor Raa (right) of 
heavy flavor electrons for Au+Au collisions at RHIC with an impact parameter of 
h = 8.2 fm together with data The cross section of gQ — > gQ is multiplied with 
the factor K = A. 



heavy flavor electron data from Ref. j4|. In BAMPS the hadronization of 
charm (bottoni) quarks to D (B) mesons is carried out with Peterson frag- 
mentation [lol . [l^ and their decay to electrons is performed by P ythia [l3l ] . 
Since the binary pQCD cross sections are too small even with running cou- 
pling and improved Debye screening [10, 14], those cross sections have been 
multiplied with a factor = 4 to be compatible with the data. Remark- 
ably, both f 2 and Raa of heavy quarks can be simultaneously described 
with K = 4. It has, of course, to be shown that currently missing contri- 
butions such as radiative processes can account for this additional factor. 
In Sec. m we will present a first attempt in this direction by considering the 
additional process gQ gQg- 

At LHC for the first time it is possible to reconstruct D mesons and, 
therefore, distinguish between charm and bottom quarks. In Fig. [2] our 
results on D mesons (l5l | is compared to data from ALICE. For the same 
parameters, that describe the RHIC data, a good agreement is also found 
at LHC. In addition, muon data from charm and bottom quarks at forward 
rapidity is also well described, as is shown in the left hand side of Fig. [3l 



4. Radiative processes 



We calculated the Gunion-Bertsch [17| matrix element for heavy quark 
scattering with a light quark within QCD to be fl8| 



\-MqQ^qQg\ = I2g^ |-^(?Q^</q| 



+ 



q± 



k\ + x2M2 (q^ - kx) 
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Figure 2. Preliminary BAMPS results on elliptic flow V2 (left) and nuclear modifi- 
cation factor Raa (right) of D mesons at Pb+Pb collisions at LHC with an impact 
parameter b together with data [1^ . The cross section of gQ — s- gQ is multiplied 
with the factor if = 4. 




Figures. Preliminary results obtained with BAMPS. Left: Raa of muons at LHC 
with data For heavy quarks only binary collisions are switched on, which are 
multiplied with K — 4. Right: Raa of heavy flavor electrons at RHIC (cf. r.h.s 
of Fig. [T]). In addition to binary collisions, radiative processes are included. No K 
factor is employed. 



This is in accordance with the result obtained with scalar QCD [19] as well 
as the result in Feynman gauge @] in the appropriate limit. Within the 
employed approximations (cf. Ref. jT3]) the matrix element of the process 
gQ gQg only differs by a color factor. 

In the right panel of Fig. [3] the Raa of heavy flavor electrons at RHIC 
is shown, which is obtained with both binary and radiative processes for 
heavy quarks without any artificial K factor. The curve is in good agree- 
ment with the data, in particular for large transverse momentum pj- where 
Peterson fragmentation works best. For smaller px coalescence might play 
an important role. 
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5. Summary 

A full space-time simulation of the quark gluon plasma phase in heavy 
ion collision is presented. The elliptic flow and nuclear modification factor of 
heavy quarks cannot be explained by binary collisions which are calculated 
with pQCD, even if a running coupling and improved Debye screening is 
employed. To be compatible with experimental data at RHIC and LHC a 
K = 4 factor is needed, which reflects that higher order contributions are 
currently missing in our model. First preliminary calculations with radiative 
contributions to heavy quark scatterings, however, look very promising since 
no artificial K factor is needed anymore to be compatible with Raa data at 
RHIC. In the future, we will investigate the impact of the implementation 
of the radiative processes on heavy flavor phenomenology in more detail. 

Acknowledgements 

J.U. would like to thank P. Gossiaux and J. Aichelin for insightful discus- 
sions. The BAMPS simulations were performed at the Center for Scientific 
Computing of the Goethe University Frankfurt. This work was supported 
by the Helmholtz International Center for FAIR within the framework of 
the LOEWE program launched by the State of Hesse. 

References 

[1] J. Adams et al. (STAR), Nucl. Phys. A757, 102 (2005), 
|5ilcl-ex/0501 009; K. Adcox et al. (PHENIX), Nucl. Phys. A757, 184 
(2005), |nucl-ex/0410003, 

[2] J. Uphoff, O. Fochler, Z. Xu, and C. Greiner, Phys. Rev. C82, 044906 
f2010). 11003 .42001 

[3] B. I. Ab elev et al. (STAR), Phys. Rev. Lett. 98, 192301 (2007), 
|iiucl-ex/0 607012; A. Adare et al. (PHENIX), Phys. Rev. Lett. 98, 
172301 (2007), ,nucl-ex/0611018l 

[4] A. Adare et al. (PHENIX), Phys.Rev. C84, 044905 ( 2011). {TUUETTETn 

[5] S. Masciocchi (ALICE), J.Phys.G G38, 124069 f2011). ITT09 . 64361 

[6] A. Rossi (ALICE), J.Phys.G G38, 124139 (2011), 1 106.59311 

[7] Y. L. Dokshitzer and D. E. Kharzeev, Phys. Lett. B519, 199 (2001), 
|hep-ph/0106202 

[8] R. Abir, C. Greiner, M. Martinez, M. G. Mustafa, and J. Uphoff (2011), 
ITi09. 55391 



6 



REFERENCES 



[9] N. Armesto, M. Cacciari, A. Dainese, C. A. Salgado, and U. A. 
Wiedemann, Phys. Lett. B637, 362 (2006), liep-ph/0511257; H. van 
Hees, V. Greco, and R. Rapp, Phys. Rev. C73, 034913 (2006), 
nucl-th/d5d8055; G. D. Moore and D. Teaney, Phys. Rev. C71, 
064904 (2005), hep-ph/0412346; M. G. Mustafa, Phys. Rev. 
C72, 014905 (2005), |he p-ph/0412402| S. Wicks, W. Horowitz, 
M. Djordjevic, and M. Gyulassy, Nucl. Phys. A784, 426 (2007), 
^cl-th/0512076; A. Adil and I. Vitev, Phys. Lett. B649, 139 (2007), 
|hep-ph/0611109; S. Peigne and A. Peshier, Phys. Rev. D77, 114017 
(2008), 0802.4364; P. B. Gossiaux and J. Aichehn, Phys. Rev. C78, 
014904 (2008), 0802.2525; W. Alberico, A. Beraudo, A. De Pace, 

A. Molinari, M. Monteno, et al., Eur.Phys.J. C71, 1666 (2011), 
11101.6 008: S. Cao and S. A. Bass (2011), 11108.51011 C. Young, 

B. Schenke, S. Jeon, and C. Gale (2011), 11111.06471 

[10] J. Uphoff, O. Fochler, Z. Xu, and C. Greiner, Phys.Rev. C84, 024908 
(2011),ill04^2295j 

[11] Z. Xu and C. Greiner, Phys. Rev. C71, 064901 (2005), 
hep-ph/0406278; Z. Xu and C. Greiner, Phys. Rev. C76, 024911 
(2007), hep-ph/0703233, 

[12] C. Peterson, D. Schlatter, I. Schmitt, and P. M. Zerwas, Phys. Rev. 
D27, 105 (1983). 

[13] T. Sjostrand, S. Mrenna, and P. Z. Skands, Comput. Phys. Commun. 
178, 852 (2008). 107 10. 38201 

[14] J. Uphoff, O. Fochler, Z. Xu, and C. Greiner, J. Phys. Conf. Ser. 
230, 012004 (2010), 1004.4091; J. Uphoff, O. Fochler, Z. Xu, and 

C. Greiner, J. Phys. Conf. Ser. 270, 012028 (2010), 1008. 1995; J. Up- 
hoff, O. Fochler, Z. Xu, and C. Greiner, Nucl.Phys. A855, 444 (2011), 
[Toil. 6183; J. Uphoff, K. Zhou, O. Fochler, Z. Xu, and C. Greiner, 
PoS BORMIO2011, 032 f201li. .1104. 2437t 

[15] O. Fochler, J. Uphoff, Z. Xu, and C. Greiner, J.Phys.G G38, 124152 
(2011), 110 7.0130 , 

[16] C. Bianchin (ALICE) (2011), these proceedings. 11111. 688B1 

[17] J. Gunion and G. Bertsch, Phys.Rev. D25, 746 (1982). 

[18] J. Uphoff, O. Fochler, Z. Xu, and C. Greiner (2011), in preparation. 

[19] P. Gossiaux, J. Aichelin, T. Gousset, and V. Guiho, J.Phys.G G37, 
094019 f2010l. [TOOT. 41661 



